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<^ ■ We present a comparison of barstrength Qt, and circumnuclear dust morphol- 

ep • ogy for 75 galaxies in order to investigate how bars affect the centers of galaxies. 

^ ■ We trace the circumnuclear dust morphology and amount of dust structure with 

structure maps generated from visible-wavelength HST data, finding that tightly 
wound nuclear dust spirals are primarily found in weakly barred galaxies. While 
strongly barred galaxies sometimes exhibit grand design structure within the cen- 
q ■ tral 10% of D25, this structure rarely extends to within ~ 10 pc of the galaxy 

nucleus. In some galaxies, these spiral arms terminate at a circumnuclear star- 
burst ring. Galaxies with circumnuclear rings are generally more strongly barred 
Q H - than galaxies lacking rings. Within these rings, the dust structure is fairly smooth 

and usually in the form of a loosely wound spiral. These data demonstrate that 
■ multiple nuclear morphologies are possible in the most strongly barred galax- 

ies: chaotic central dust structure inconsistent with a coherent nuclear spiral, a 
grand design spiral that loses coherence before reaching the nucleus, or a grand 
^ ■ design spiral that ends in a circumnuclear ring. These observations may indicate 

that not all strong bars are equally efficient at fueling material to the centers 
of their host galaxies. Finally, we investigate the longstanding hypothesis that 
SB(s) galaxies have weak bars and SB(r) galaxies have strong bars, finding the 
opposite to be the case: namely, SB(r) galaxies are less strongly barred and have 
less dust structure than SB(s) galaxies. In general, more strongly barred galaxies 
tend to have higher nuclear dust contrast. 

Subject headings: dust - galaxies: ISM - galaxies: nuclei - galaxies: spiral - 
galaxies: structure - ISM: structure 



1 Based on observations with the NASA/ESA Hubble Space Telescope obtained at the the Space Tele- 
scope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Incorporated, under NASA contract NAS5-26555. 
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1. Introduction 

Most spiral galaxies have large-scale bars, and any complete picture of galaxies must 
include the impact of bars on their evolution. Bars are important because they are an obvious 
mechanism for funneling gas and dust to the centers of galaxies: gas clouds orbiting in the 
disk lose angular momentum as they encounter the bar, thereby sinking towards the galaxy's 
center (Binney & Tremaine 1987; Athanassoula 1992; Piner et al. 1995; Regan et al. 1997; 
Maciejewski et al. 2002). Bars are therefore expected to play a major role in circumnuclear 
star formation, bulge growth, and the fueling of the central, supermassive black hole. 

Bar-driven radial gas inflow should drive large quantities of gas and dust into the central 
regions of galaxies, thereby increasing the central gas concentration of barred galaxies relative 
to unbarred galaxies. Observations by Sakamoto et al. (1999) and Sheth et al. (2005) support 
this view, finding higher concentrations of molecular gas in the central kiloparsec of barred 
galaxies than in unbarred galaxies. As higher gas density is empirically correlated with 
higher star formation rates (e.g., Kennicutt 1998), there should be a correlation between 
bar-driven inflow and central star formation. Observational evidence for enhanced nuclear 
or circumnuclear star formation in barred galaxies provides broad support for this picture 
(Ho et al. 1997; Maoz et al. 2001; Knapen et al. 2006). 

It has long been speculated that this same radial gas inflow could provide active galactic 
nuclei (AGN) with fuel (Simkin et al. 1980; Schwarz 1981). There is so far no evidence, 
though, that bars (large-scale or nuclear) are the primary fuel source for AGN: not all active 
galaxies are barred, not all barred galaxies are active, and, in fact, well-matched samples of 
active and inactive galaxies have the same bar fraction (Ho et al. 1997; Mulchaey & Regan 
1997; Laine et al. 2002). Ho et al. speculate that this is because radially-transported gas is 
prevented from reaching the nucleus. 

If bars are indeed fueling galaxies' circumnuclear regions, then this should be reflected in 
the morphology of the cold circumnuclear interstellar medium (ISM). Martini et al. (2003a) 
establish a nuclear dust morphology classification system with four spiral classes (grand 
design, tightly wound, loosely wound, and chaotic spiral) and two non-spiral classes (chaotic 
and no structure) to quantify differences between various galaxy types. Martini et al. (2003b) 
find no differences in the nuclear dust structure of active and inactive galaxies. This study 
does find evidence that nuclear grand design spirals are primarily found in barred galaxies 
and that these grand design spirals connect to large scale dust lanes in barred galaxies. 
Likewise, Martini et al. primarily find tightly wound nuclear spirals in unbarred galaxies. 

All of these studies of the role of bars in fueling central star formation or black holes sim- 
ply compare "barred" versus "unbarred" galaxies. This discretization glosses over, and over 
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simplifies, the long- known continuum of barstrengths (de Vaucouleurs 1959). One early way 
of quantifying the strength of a bar is the deprojected bar ellipticity e b (Martin 1995). The 
disadvantage of using the bar ellipticity is that defining the bar — i.e., its exact dimensions — is 
a somewhat subjective process (Buta & Block 2001). Abraham & Merrifield (2000) introduce 
the parameter /bar, which is a function of the bar axis ratio, and therefore susceptible to the 
same systematics as e b . Buta & Block, by expanding on the method of Combes & Sanders 
(1981), circumvent this problem by using a force ratio, which they call Q b . Assuming that 
the light traces the underlying mass distribution, they use near-infrared images to calcu- 
late a map of a galaxy's potential. From these potential maps, and with some assumptions 
about the characteristic scale-height, they then calculate the tangential force F T and mean 
(axisymmetric) radial force (Fr). They define the ratio map Qt as 



This ratio map has the property that in each of four quadrants Qt reaches a local ex- 
tremum. Letting Q^ ax,k be the absolute value of such an extremum in the kth quadrant, 
the barstrength Q b is defined as 



The cited uncertainties on Q b are a measure of how much Q™ ax varies by quadrant. From 
a comparison with the three barstrength classes defined by de Vaucouleurs (unbarred SA, 
weakly barred SAB, and strongly barred SB), Buta & Block find a reasonable correlation 
between the RC3 bar classification and Q b , although with substantial scatter. In particular, 
SA galaxies have Q b < 0.1, SAB galaxies have 0.05 < Q b < 0.2, and SB galaxies have 
Q b > 0.15. The most strongly barred galaxies have Q b ~ 0.6 (Buta & Block 2001). While 
we adopt Q b as the best available measure of barstrength, we note that Athanassoula (1992) 
shows that both the quadrupole moment (strength) and the pattern speed play an important 
role in how effectively the bar drives mass towards the center of the galaxy. Pattern speed, 
however, is much more difficult to measure than parameters that can be calculated from 
photometric data. 

Martini (2004) used Q b and the nuclear classifications of Martini et al. (2003a) to com- 
pare the barstrengths and nuclear dust morphologies of 48 galaxies with archival Hubble 
Space Telescope (HST) data. He found that grand design nuclear spirals are primarily found 
in strongly barred galaxies, while axisymmetric tightly wound nuclear spirals are exclusively 
found in galaxies with Q b < 0.1. With a larger sample and reconsideration of the circum- 
nuclear dust morphology classification at the smallest scales, we examine how nuclear dust 
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structure varies with barstrength. We identify 75 galaxies, described in §2, with a measured 
barstrength Qf, from the literature. To analyze the dust morphology, we create "structure 
maps" from archival HST data using the technique of Pogge & Martini (2002) as discussed in 
§3.1. We then classify the galaxies according to a refined version of the classification scheme 
proposed in Martini et al. (2003a), as described in §3.2. The main difference between the 
classification scheme used here and that of Martini et al. (2003a) is a stronger focus on the 
central-most regions of the galaxy. We discuss our results in §4. 



2. Data 

We began our study with a data set of 127 galaxies with both a measured barstrength 
(primarily from the Ohio State University Bright Galaxy Survey (OSUBGS) by Laurikainen 
et al. 2004) and archival HST data. The HST images were obtained with WFPC2, with the 
exception of three images from the Advanced Camera for Surveys (ACS), and most employed 
either the F606W or F814W filters. The galaxies studied by the OSUBGS were selected to 
be a representative sample of the nearby universe; specifically, based on the RC3 T-type 
distribution (0 < T < 9, or SO/a to Sm) to ensure there is no morphological selection bias 
(Eskridge et al. 2002). 
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Table 1. Data 



Galaxy 


Nuclear 




Qb 




o"sm 


Detector 


Filter 


HST 


T 


D 


r c 


Notes 




Class 








(xlO- 


) 






Prop. ID 




(Mpc) 


(pc) 




(1) 


(2) 




(3) 






V 1 ) 




(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC289 


LW 


0.212 


± 


0.003 


1.280 


± 


0.030 


PCI 


F606W 


6359 


4 


19.4 


221 d 


SB(rs) 


NGC488 


LW 


0.032 


± 


0.003 


0.951 


± 


0.015 


PCI 


F606W 


6359 


3 


29.3 


447 




NGC613 


cs 


0.401 


± 


0.045 


4.469 


± 


0.017 


WF3 


F606W 


9042 


4 


17.5 


280 


SB(rs) 


NGC864 


cs 


0.360 


± 


0.037 


3.726 


± 


0.064 


PCI 


F606W 


8597 


5 


20 


177 


LGD 


NGC972 


c 


0.220 


± 


0.060 a 


3.948 


± 


0.060 


PCI 


F606W 


6359 


2 


21.4 


206 




NGC1068 


c 


0.165 


± 


0.010 


4.204 


± 


0.072 


PCI 


F606W 


8597 


3 


14.4 


134 d 




NGC1073 


c 


0.607 


± 


0.013 


1.839 


± 


0.013 


WF3 


F606W 


9042 


5 


15.2 


217 


SB(rs) 


NGC1097 


LW 


0.279 


± 


0.048 


2.416 


± 


0.025 


HRC 


F814W 


9788 


3 


14.5 


394 


SB(s), LGD 


NGC1241 


LW 


0.251 


± 


0.028 


3.383 


± 


0.037 


PCI 


F606W 


5479 


3 


46.6 


382 


SB(rs), LGD 


NGC1300 


LW 


0.537 


± 


0.011 


2.362 


± 


0.025 


PCI 


F606W 


8597 


4 


18.8 


311 d 


SB(rs), LGD 


NGC1317 


CS 


0.085 


± 


0.007 


1.910 


± 


0.023 


PCI 


F606W 


5446 


1 


16.9 


135 


LGD 


NGC1326 


cs 


0.160 


± 


0.020 b 


2.506 


± 


0.019 


PCI 


F555W 


6496 


-1 


16.9 


191 


LGD 


NGC1350 


c 


0.243 


± 


0.039 


0.863 


± 


0.008 


PCI 


F606W 


5446 


1.8 


16.9 


258 


SB(r) 


NGC1365 


cs 


0.400 


± 


0.110 a 


3.892 


± 


0.008 


PCI 


F606W 


8597 


3 


16.9 


552 


SB(s), LGD 


NGC1385 


c 


0.319 


± 


0.030 


2.298 


± 


0.024 


PCI 


F606W 


6359 


6 


17.5 


172 


SB(s) 


NGC1398 


N 


0.202 


± 


0.011 


0.552 


± 


0.001 


PCI 


F606W 


8597 


2 


16.1 


332 


SB(r) 


NGC1433 


GD 


0.370 


± 


0.060 a 


1.772 


± 


0.028 


WF3 


F814W 


9042 


2 


11.6 


218 


SB(r) 


NGC1530 


TW 


0.610 


± 


0.160 a 


3.965 


± 


0.005 


PCI 


F606W 


8597 


3 


36.6 


487 


SB(rs), LGD 


NGC1637 


C 


0.202 


± 


0.014 


2.054 


± 


0.061 


WF4 


F555W 


9155 


5 


8.9 


91 d 


LGD 


NGC1808 


LW 


0.274 


± 


0.001 


5.326 


± 


0.046 


PCI 


F675W 


6872 


1 


10.8 


203 




NGC2196 


C 


0.070 


± 


0.005 


0.858 


± 


0.017 


PCI 


F606W 


6359 


1 


28.8 


236 




NGC2442 


CS 


0.669 


± 


0.428 


2.559 


± 


0.031 


WF3 


F606W 


9042 


3.7 


17.1 


273 




NGC2655 


cs 


0.128 


± 


0.004 


1.822 


± 


0.003 


PCI 


F547M 


5419 





24.4 


348 




NGC2775 


N 


0.050 


± 


0.010 


0.579 


± 


0.009 


PCI 


F606W 


6359 


2 


17 


211 




NGC2964 


c 


0.310 


± 


0.003 


6.463 


± 


0.162 


PCI 


F606W 


6359 


4 


21.9 


149 d 


LGD 


NGC2985 


TW 


0.056 


± 


0.001 


1.608 


± 


0.033 


WF2 


F606W 


5479 


2 


22.4 


298 




NGC2997 


TW 


0.060 


± 


0.020 b 


3.498 


± 


0.067 


HRC 


F555W 


9989 


5 


13.8 


358 




NGC3054 


C 


0.170 


± 


0.020 b 


1.697 


± 


0.013 


PCI 


F606W 


6359 


3 


26.5 


207 




NGC3077 


C 


0.119 


± 


0.016 


6.253 


± 


0.093 


WF3 


F814W 


9144 


90 


2.1 


33 




NGC3081 


GD 


0.170 


± 


0.020 b 


2.501 


± 


0.023 


PCI 


F606W 


5479 





32.5 


198 


LGD 


NGC3169 


C 


0.090 


± 


0.005 


3.254 


± 


0.066 


PCI 


F547M 


5149 


1 


19.7 


250 




NGC3227 


C 


0.158 


± 


0.021 


4.615 


± 


0.109 


PCI 


F606W 


5479 


1 


20.6 


181 d 




NGC3310 


C 


0.060 


± 


0.010 c 


1.575 


± 


0.007 


PCI 


F814W 


6639 


4 


18.7 


168 




NGC3338 


TW 


0.083 


± 


0.005 


1.318 


± 


0.011 


WF3 


F606W 


9042 


5 


22.8 


391 




NGC3359 


C 


0.460 


± 


0.050 c 


1.914 


± 


0.002 


WF3 


F606W 


9042 


5 


19.2 


405 


SB(rs) 


NGC3486 


GD 


0.108 


± 


0.002 


0.968 


± 


0.006 


PCI 


F606W 


8597 


5 


12.3 


165 d 




NGC3504 


LW 


0.288 


± 


0.030 


3.862 


± 


0.125 


PCI 


F606W 


5479 


2 


26.5 


207 




NGC3898 


C 


0.047 


± 


0.000 


0.877 


± 


0.009 


PCI 


F606W 


6359 


2 


21.9 


278 




NGC4027 


C 


0.623 


± 


0.008 


2.143 


± 


0.021 


PCI 


F814W 


8599 


8 


25.6 


175 d 


SB(s) 


NGC4030 


TW 


0.060 


± 


0.013 


2.193 


± 


0.026 


PCI 


F606W 


6359 


4 


25.9 


314 




NGC4051 


CS 


0.280 


± 


0.008 


2.873 


± 


0.082 


PCI 


F606W 


5479 


4 


17 


260 




NGC4138 


c 


0.046 


± 


0.007 


2.829 


± 


0.008 


PCI 


F547M 


6837 


-1 


17 


127 




NGC4254 


LW 


0.122 


± 


0.029 


2.640 


± 


0.029 


PCI 


F606W 


8597 


5 


16.8 


262 




NGC4303 


LW 


0.259 


± 


0.044 


3.928 


± 


0.062 


HRC 


F555W 


9776 


4 


15.2 


218 d 


LGD 
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Table 1 — Continued 



Galaxy 


Nuclear 




Wb 




usm 


Detector 


r liter 


no 1 


T 
1 


LJ 


Tc 


IN OtCS 












(xlO- 








riop. 11J 




(ivipcj 


(P C ) 






( z ) 










(4) 






W 


(I) 


\°) 






Cl 1 "i 


NGC4314 


LW 


0.442 


± 


0.024 


1.390 


± 


0.029 


PCI 


F606W 


8597 


i 


9.7 


118 


SB(rs), LGD 


NGC4321 


LW 


0.183 


± 


0.027 


2.123 


± 


0.045 


PCI 


F702W 


5195 


4 


16.8 


362 


LGD 


NGC4414 


LW 


0.149 


± 


0.003 


1.857 


± 


0.032 


PCI 


F606W 


8597 


5 


9.7 


118 




NGC4450 


C 


0.131 


± 


0.011 


0.634 


± 


0.005 


PCI 


F814W 


5375 


2 


16.8 


256 




NGC4501 


LW 


0.072 


± 


0.026 


1.834 


± 


0.023 


PCI 


F606W 


6359 


3 


16.8 


276 d 




NGC4504 


C 


0.136 


± 


0.019 


1.723 


± 


0.016 


WF3 


F814W 


9042 


6 


19.5 


248 




NGC4579 


cs 


0.197 


± 


0.020 


1.713 


± 


0.032 


PCI 


F547M 


6436 


3 


16.8 


288 




NGC4593 


LW 


0.309 


± 


0.020 


1.963 


± 


0.069 


PCI 


F606W 


5479 


3 


39.5 


447 


SB(rs), LGD 


NGC4651 


GD 


0.120 


± 


0.046 


1.329 


± 


0.022 


PCI 


F555W 


5375 


5 


16.8 


195 




NGC4698 


N 


0.084 


± 


0.040 


0.700 


± 


0.006 


PCI 


F606W 


6359 


2 


16.8 


195 




NGC4736 


LW 


0.048 


± 


0.004 


1.684 


± 


0.003 


PCI 


F555W 


5741 


2 


4.3 


140 




NGC4939 


C 


0.128 


± 


0.052 


2.297 


± 


0.008 


PCI 


F606W 


5479 


4 


17.3 


113 




NGC4941 


CS 


0.056 


± 


0.008 


1.582 


± 


0.017 


PCI 


F606W 


8597 


2 


44.3 


531 d 


LGD 


NGC5054 


cs 


0.090 


± 


0.023 


1.790 


± 


0.032 


PCI 


F606W 


8597 


4 


9.7 


102 




NGC5121 


GD 


0.024 


± 


0.007 


3.457 


± 


0.003 


PCI 


F606W 


6359 


1 


27.3 


340 d 




NGC5194 


LW 


0.160 


± 


0.000 c 


1.198 


± 


0.006 


PCI 


F547M 


5123 


4 


22.1 


125 




NGC5236 


C 


0.190 


± 


0.040 b 


1.942 


± 


0.027 


PCI 


F814W 


8234 


5 


7.7 


251 




NGC5248 


TW 


0.269 


± 


0.064 


1.717 


± 


0.027 


PCI 


F814W 


6738 


4 


2.7 


92 d 




NGC5427 


LW 


0.231 


± 


0.074 


2.327 


± 


0.010 


PCI 


F606W 


5479 


5 


22.7 


376 d 


LGD 


NGC5643 


GD 


0.415 


± 


0.013 


1.150 


± 


0.007 


PCI 


F606W 


8597 


5 


38.1 


312 


LGD 


NGC6217 


C 


0.360 


± 


0.010 c 


6.472 


± 


0.171 


PCI 


F606W 


5479 


4 


23.9 


210 


SB(rs) 


NGC6221 


C 


0.436 


± 


0.112 


5.053 


± 


0.010 


PCI 


F606W 


5479 


5 


19.4 


200 


SB(s) 


NGC6300 


CS 


0.187 


± 


0.002 


2.505 


± 


0.074 


PCI 


F606W 


5479 


3 


14.3 


186 


SB(rs), LGD 


NGC6384 


c 


0.136 


± 


0.020 


0.910 


± 


0.009 


PCI 


F606W 


6359 


4 


26.6 


477 




NGC6814 


GD 


0.070 


± 


0.010 a 


1.585 


± 


0.049 


PCI 


F606W 


5479 


4 


22.8 


205 




NGC6951 


LW 


0.280 


± 


0.040 a 


2.036 


± 


0.041 


PCI 


F606W 


8597 


4 


24.1 


273 


LGD 


NGC7098 


C 


0.200 


± 


0.020 b 


0.856 


± 


0.003 


PCI 


F555W 


6633 


1 


29.1 


345 




NGC7213 


TW 


0.023 


± 


0.002 


1.124 


± 


0.020 


PCI 


F606W 


5479 


1 


22 


198 




NGC7479 


C 


0.696 


± 


0.060 


2.022 


± 


0.046 


PCI 


F814W 


6266 


5 


32.4 


384 


SB(s) 


NGC7552 


LW 


0.395 


± 


0.044 


5.391 


± 


0.267 


PCI 


F606W 


5479 


2 


19.5 


192 


SB(s) 


NGC7727 


LW 


0.096 


± 


0.024 


1.222 


± 


0.031 


PCI 


F555W 


7468 


1 


23.3 


317 





Note. — iterature data and measurements for the 75 galaxies in our sample. The classification codes in column 2 are: GD: 
grand design nuclear spiral; TW: tightly wound spiral; LW: loosely wound spiral; CS: chaotic spiral; C: chaotic circumnuclcar 
dust; N: no circumnuclear dust structure. These codes are defined in §3.2. Qi, (col. [3]) is from Laurikainen et al. (2004) unless 
otherwise noted. See §3.3 for a discussion of the structure map rms <x sm (col. [4]). In column 5, "HRC" is the High Resolution 
Channel for ACS; the other detectors refer to different chips on WFPC2. The Hubble T-types (col. [8]) and SB classifications 
(col. [11]) arc from the RC3. The distances (col. [9]) are from Tully & Fisher (1988) and assume H = 75 km s _1 Mpc -1 . 
The solid white circles in Fig. 1 have radius r c (col. [10]) as defined in eqn. (3). LGD (col. [11]) refers to "large grand design" 
structure; see §4.1 for details. 

a Block ct al. (2004) 

b Buta & Block (2001) 

c Laurikainen & Salo (2002) 
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d r c = r eff (bulge) 
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For each galaxy, we define a critical "central" radius r c by 

r c = min{r cff (bulge), 0.01 • D 25 }, (3) 

where r e fj(bulge) is the bulge radius reported in Laurikainen et al. (2004) and D 2 s is the 
RC3 value of the 25th magnitude B-band isophotal diameter (de Vaucouleurs et al. 1991). 
If Laurikainen et al. give no bulge radius, or if the galaxy is completely bulgeless, then 
r c = 0.01-D 2 5- As the bulge will dominate the circumnuclear region, the minimum is used in 
the definition of r c to guarantee that the central region is no larger than the bulge. Most 
(about 80%) of the galaxies in the final sample have r c = 0.01-D 2 5- A circle of radius r c , which 
corresponds to 30-550 pc projected, is included on the structure maps used for classification 
(see §3.2). 

Objects with high inclination (axis ratio i? 2 5 < 0.30) or low signal-to-noise ratio (< 10) 
are excluded from the sample. We also discard low-resolution objects, given by r c < 20 res- 
olution elements, as well as ones with unfavorable location on the chip or highly saturated 
centers. In general, the objects failing the resolution or signal-to-noise ratio cuts have no 
clearly discernible coherent spiral arm structure. All objects are closer than 50 Mpc. Our 
final sample, summarized in Table 1, consists of the 75 galaxies which pass our resolution and 
signal-to-noise ratio cuts; the T-type distribution of the Eskridge et al. sample is roughly 
maintained. 



3. Analysis 

We employ two methods of examining the nuclear dust content in galaxies. The first 
is to look at what structures the dust forms, i.e., a morphological classification, and the 
second is to simply study the amount of dust structure in the central regions. We describe 
in §3.1 how we use the structure map technique of Pogge & Martini (2002) to enhance 
dust structures on the scale of 1-15 pc. We then use these structure maps to classify the 
circumnuclear dust morphology into six classes (four types of spirals, plus chaotic and no 
structure), as described in §3.2. In §3.3 we describe how we estimate the relative amounts of 
central dust structure in different galaxies. Finally, in §3.4, we discuss how we statistically 
compare pairs of distributions of galaxies. 

3.1. Structure Maps 

As it is difficult to obtain a large sample of high quality, high resolution multi-band 
(e.g., both V and H) data, we used the "structure map" technique developed by Pogge & 
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Martini (2002) to enhance the observable dust structure in each galaxy. Structure maps 
can be thought of as "single-band color maps;" while requiring data in only one band, 
they provide a sharp contrast between dust and star-forming regions. Martini & Pogge 
(1999) note that the regions of highest contrast in V — H color maps are primarily due to 
absorption by dust structures in the V image alone, while the smoother if -band images trace 
the underlying stellar population. Therefore, the i/-band images serve mostly to subtract 
larger-scale surface brightness variations from the ^-band images in the V — H color maps. 
By making structure maps from images taken at visible wavelengths, it is possible to enhance 
the structure due to dust and emission regions. Furthermore, as visible-wavelength images 
from HST have higher resolution than those in the near infrared (both being diffraction 
limited), a structure map can make the most of high quality data. 

The structure map is mathematical defined as 

where S is the structure map, I is the image, P is the point spread function (PSF), P l 
is the transform of the PSF, and <g> is the convolution operator. Structure maps bring out 
variations on the smallest resolvable scale of the image, i.e., that of the PSF. For our objects, 
this corresponds to a projected scale of l-15pc. The PSFs are modelled using the TinyTim 
software (Krist & Hook 1999). 

Formally, structure maps are similar to the second-order iteration of the Richardson- 
Lucy (R-L) image correction (Richardson 1972; Lucy 1974). At this step in the image 
restoration, the first-order, smooth structure of the image has been removed, but the struc- 
tures on the scale of the size of the PSF remain. Structure maps for the 75 galaxies in our 
sample are given in Figure 1; in all of the structure maps presented in this paper, the dusty 
regions appear dark while enhanced stellar light and emission-line regions are bright (e.g., 
the F606W filter admits several bright emission lines.) 
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Fig. 1. — Structure maps for 75 galaxies with measured barstrengths. Each panel shows 
the inner 5% of -D25 from the RC3 catalog. The solid white circles are given by a radius of 
r c , as defined in eq. (3). In the case that r c g(bulge) > r c , the bulge radius is marked by a 
dashed white circle; these dashed circles were not included on the structure maps used for 
classification, but are merely given here for reference. Dark regions are due to dust, while 
bright regions are due to emission. Several images (e.g., NGC2997) also include the ACS 
coronographic finger. North is up and East is to the left. 




Fig. 1. — Continued. 




Fig. 1. — Continued. 
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3.2. Nuclear Dust Classification 

We classify the galaxies according to a refined version of the nuclear classification scheme 
presented in Martini et al. (2003a). The main distinction between the scheme used here and 
the original one is that we have a stronger emphasis on the central-most regions of the galaxy; 
other differences between the scheme used here and the one in Martini et al. (2003a) are 
discussed below. The structure maps used for classification are scaled to have a width and 
height equal to 5% of D 2 5- The classification scheme uses six categories: 

Grand Design (GD): Two spiral arms symmetric about a 180° rotation, at least one of 
which extends all the way to the unresolved nucleus of the galaxy, and at least one of 
which is a dominant feature 

Tightly Wound (TW): Coherent spiral arm structure over a large range in radius; pitch 
angle less than or equal to 10° 

Loosely Wound (LW): Coherent spiral arm structure over a large range in radius; pitch 
angle greater than 10° 

Chaotic Spiral (CS): Unambiguous evidence for spiral arm structure with a unique sense 
of chirality, but not coherent over a large range in radius 

Chaotic (C): Dust structure with no evidence of spiral structure 

No Structure (N): No discernible dust structure within the central region 

Examples of each classification are given in Figure 2, and the classifications for all of 
the galaxies in our sample are given in Table 1. There are several subtle, yet important, 
refinements to the original classification scheme. First, the classification is based on a region 
that scales with the size of the galaxy (the central 5% of D25) rather than within the central 
19.6" used by Martini et al. (2003a). Most changes in classification can be attributed to this 
more physically defined focus. The largest difference in an individual class is the requirement 
that for an object to be classified as GD, at least one of the arms must be observed to extend 
to the nucleus of the galaxy. For this reason, some galaxies previously classified as GD now 
have a different classification (e.g., NGC1365). We comment on this distinction in greater 
detail in §4.1. 



-14- 




Fig. 2. — Structure maps of prototypes for the six circumnuclear morphology classes dis- 
cussed in §3.2. The white circles are given by a radius of r c , as defined in eq. (3). Dark 
regions are dust, while bright regions are emission. 
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We also quantify the distinction between tightly and loosely wound circumnuclear spi- 
rals by setting 10° as the delineation between "small" and "large" pitch angles. This division 
roughly corresponds to that of Sb galaxies (Kennicutt 1981). Previously, the differentiation 
between TW and LW was done purely by eye, and also took into account arm morphology. 
Specifically, for a nuclear spiral to be classified as "tightly wound," the arms had to be 
traceable through at least one complete revolution, and therefore the classification was po- 
tentially subject to signal-to-noise ratio bias (for an arm to remain coherent about an entire 
revolution, and thus be classified as TW, a higher SNR was needed). We measure pitch 
angles with an interactive PGPlot program which fits logarithmic spirals to a deprojected 
image of the structure map. Five points of increasing radius on a single prominent spiral 
arm of the galaxy are chosen within an annulus given by 25-75% of the bulge radius, or of 
1% of D 2 5 if the galaxy is either bulgeless or there is no available effective bulge radius data. 1 
The program then fits logarithmic spirals to each consecutive pair of points, averaging the 
resulting four pitch angles to give the pitch angle for the entire arm. As found by Kennicutt 
(1981) for large-scale spirals, the nuclear spirals are not well fit by logarithmic spirals — a 
given arm does not have constant pitch angle as a function of radius. Still, this method 
yields an "average" pitch angle for the arm. The observed non-constant pitch angle along a 
spiral arm implies that the region is experiencing differential rotation which is inconsistent 
with constant rotational velocity (Maciejewski 2004). We also note that the pitch angle often 
varies between different arms in a single galaxy. Notwithstanding these caveats, we find that 
this method suffices for determining whether or not the pitch angle is above or below 10°, 
thus differentiating between TW and LW nuclear spirals. 

Finally, the "central region" used in classifying the N class of objects was previously 
defined as a "few hundred parsecs" (Martini et al. 2003a). We have clarified this so that 
the central region is specific to each galaxy, namely, a circle of radius r c . While in general 
the entire region within 0.05D 2 5 is used for classification, the region within r c is also used 
as a guide in deciding between different classifications; if a galaxy is seen to have a different 
morphology at small radii than at large radii, such as is seen in NGC1068, the classification at 
small radii is used. This decision also led to changes in the classifications of several galaxies 
from those reported in previous papers. Though this distinction is subject to the definition 
of r c , only two galaxies, NGC1068 and NGC3227, would have their classifications changed 
(LW to C and CS to C, respectively) if r c were defined to be 1% of D25 for all galaxies, 
rather than the minimum used in Equation (3). Despite this, we choose to keep the more 



lr rhis radius is distinct from r c in that if Laurikainen et al. (2004) report a non-zero bulge radius, then 
it is used, rather than taking the minimum of r e g-(bulge) and 0.01£>25- This is because at larger radii, there 
is more likely to be a coherent arm with a measurable pitch angle. 
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physically motivated definition of r c for reasons discussed in §2. 

3.3. Dust Contrast 

To quantify the amount of relative structure galaxies, we calculated the root mean 
squared (rms) of the pixel-to-pixel variations of the galaxies' structure maps within r c . The 
more dust structure or star formation in a galaxy the higher the dust contrast in the structure 
map, and therefore the higher the rms. We assign an uncertainty to the rms by calculating 
the standard deviation of 11 measurements of the rms within a circular aperture ranging in 
radius from 0.95r c to 1.05r c in increments on 0.01r c . We further define <r sm to be the average 
of these 11 measured rms values. The variation on <r sm is found to be no more than 5% for 
any object. The measured rms values and uncertainties for each galaxy are given in Table 1. 

Because the structure map is sensitive to structures on the scale of the PSF, rather than 
some physical scale, there may be a strong distance-dependent bias in cr sm . We therefore 
tested how <j sm changes for a given galaxy upon changing the resolution using a test sample of 
nine galaxies with high enough resolutions to meet our resolution and signal-to-noise cuts (see 
§2) even after being convolved with a Gaussian with a standard deviation of 2.6 pixels, which 
corresponds to more than doubling their distance. For each object, we convolved both the 
image and the PSF with a Gaussian with a standard deviation of 1.0-2.6 pixels (in 0.1 pixel 
increments) , from which a set of 17 degraded structure maps were constructed. An example of 
this degradation for NGC4321 is shown in Figure 3. We then calculated <7 sm (and uncertainty, 
as described above) for each of these degraded structure maps. While <r sm did change for 
each object, the change was found to be comparable with the calculated uncertainties for 
any given structure map. This implies that at these scales, the dust structure is roughly 
scale-invariant, in agreement with the findings of Elmegreen et al. (2002) that the power 
spectra of dust structure becomes relatively shallow at small scales. We thus conclude that 
a sm does not contain a distance-dependent bias and employ it below to compare different 
subsamples of galaxies. 



Fig. 3. — Degraded structure maps of NGC4321 for testing the dependence of structure map 
rms <7 sm on resolution. The percents in the upper left-hand corners are the percent change in 
a sm relative to the original structure map. The numbers in the upper right-hand corners are 
the standard deviations in pixels of the Gaussians with which both the image and the PSF 
were convolved before construction of the structure map. The white circles have a radius of 
r c , as defined in eq. (3). The resolution of the original image is 0.06", or 4.9 pc, while the 
resolution of the 2.6-a Gaussian convolved image is 0.13", or 10.8 pc. This is equivalent to 
moving the galaxy from 16.8 Mpc to 37 Mpc. Relative to the measured uncertainty on the 
original <r sm (2.1%), there is no significant change in a sm with resolution. 
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3.4. Comparing Q b and a sm Distributions 

Once the galaxies are classified, we have sets of barstrengths Qb and structure map rms 
values <7 sm for each nuclear classification. We then want to answer the question of how likely 
it is that the objects of a given classification are drawn from the same parent distribution 
as other galaxies. The standard method for comparing two arbitrary distributions is the 
Kolmogorov-Smirnov (K-S) test, which calculates the maximum vertical separation Dks 
between two cumulative distributions. The larger the K-S distance Dks for two samples, 
the less likely the two distributions represent the same parent population. Dks only takes 
on discrete values, namely, multiples of 1/(N A + Nb), where Na and N B are the number 
of elements in the two samples A and B under comparison. To compare two samples, we 
first calculate the distance Dks between the two distributions. Then, for 10,000 trials, we 
compare this Dks to the distance between two samples of size Na and Nb drawn from a 
uniform distribution, yielding a probability that the two samples are in fact drawn from the 
same distribution. These probabilities do not change significantly when errors on Q b or cr sm 
are taken into account. In the event that the two samples intersect, we exclude the shared 
galaxies from both samples in order to do the comparison. 

In the few cases where both sample sizes small, we use the Wilcoxon rank-sum test 
instead of the K-S test. The Wilcoxon test compares two samples A and B of size N A and 
N B , with N A < N B , where the "true" mean values of A and B are ha and /j, b , respectively. 
Letting the sum of the ranks in sample A be w, the probability that \i A > Hb is equal to the 
probability that the sum of the ranks of Na randomly chosen elements (from a set of size 
Na + N B ) is < w, while the probability that /ia < I^b is equal to the probability that the 
sum of ranks is > Na(Na + N B + 1) — w. 

4. Results and Discussion 

Our classifications and measurements of cr sm are summarized in Table 1. There is a 
2% probability that the tightly wound (TW) nuclear dust spirals have the same underlying 
distribution of barstrength as the rest of the sample. As shown in Figure 4, TW nuclear 
spirals are found primarily — although not exclusively — in weakly barred galaxies. As the 
TW class is defined to be those dust spirals with pitch angle < 10°, this is consistent with 
the idea that a galaxy which is axisymmetric at large scales (i.e., has a small Qb) will either 
have high differential rotation or simply that the central regions of these galaxies cannot 
be highly disrupted (Maciejewski 2004). All four late-type galaxies in our sample (T > 6) 
have chaotic circumnuclear dust structure, as shown in the right-hand panel of Figure 5. No 
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such differences were found for any other nuclear class, 2 as shown in the left-hand panel of 
Figure 5 and Figure 6. As previous studies have suggested a connection between large-scale 
bars and grand design nuclear spirals, we begin our discussion of strongly barred galaxies 
below with a consideration of grand design nuclear spirals. 



2 While all three N class galaxies have a T = 2 (corresponding to a Hubble type of Sab), this is not 
significant. 
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Fig. 4. — Comparison of Qb for TW galaxies and all other galaxies. The vertical widths 
correspond to the central 68-percentile spread of the cumulative fraction and are due to the 
uncertainty in Q^. Galaxies with T W circumnuclear dust spirals are more weakly barred than 
typical galaxies, with a 2% probability of being drawn from the same parent distribution. 
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Fig. 5. — Left: No connection is seen between Hubble T-type for T < 5 (Sc) and nuclear 
classification. Right: Structure maps for the four galaxies in the sample with T > 6. The 
large scale morphologies from the RC3 catalog are: NGC3077, 10; NGC4504, Scd; NGC1385, 
Scd; NGC4027, Sdm. All four of these galaxies are classified as having chaotic circumnuclear 
dust structure. 




Fig. 6. — Nuclear classifications and barstrength Qt,. Left: Distribution of nuclear classifica- 
tion into different barstrength classes. Right: Monte Carlo average cumulative distribution 
of barstrength in each nuclear class and total sample. 
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4.1. Circumnuclear Grand Design Spirals 

Hydro dynamic simulations suggest that large-scale bars can lead to the formation of 
grand design nuclear spirals (Patsis & Athanassoula 2000; Englmaier & Shlosman 2000), 
and several observational studies support this hypothesis (Pogge & Martini 2002; Martini 
et al. 2003b). With our focus on the very central regions of nuclear spirals, we find that this 
is in fact not the case: the grand design nuclear spirals described in §3.2 are not preferentially 
found in strongly barred galaxies. We find, however, that symmetric two-arm spiral structure 
is common at larger scales in strongly barred galaxies. To differentiate between circumnuclear 
grand design spirals found at different scales, we hereafter refer to grand design structure at 
small radii, defined as GD in §3.2, as SGD structure (small grand design), and prominent 
grand design structure at larger radii (within 10% of -D25, as shown in Figure 7) as LGD 
structure (large grand design). 

Figure 7 shows four examples of galaxies for which grand design structure is evident at 
larger scales, but is not present at smaller scales. Seven galaxies (9% of the entire sample) 
have SGD structure, while twenty (26%) have LGD structure. This larger-scale GD structure 
extends to the center of the galaxy in only two (10%) of these twenty galaxies. As can be 
seen in Figure 8, LGD galaxies are preferentially more strongly barred than SGD galaxies. 
A K-S test excluding the two shared galaxies reveals a 3% probability that SGD and LGD 
barstrengths are drawn from the same parent distribution. We further find a 0.4% probability 
that LGD galaxies have the same underlying distribution of Qb as the other 59 galaxies in 
the sample; galaxies with LGD nuclear spirals are more strongly barred than other galaxies. 
It is interesting to note that the two LGD galaxies with the smallest Qb (NGC4941 and 
NGC1317) also have an obvious nuclear bar within the grand design structure (Erwin 2004 
and Greusard et al. 2000, respectively). 
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Fig. 7. — Four examples of galaxies with LGD (large grand design) but not SGD (small 
grand design) structure. Left: Width and height are given by 0.1D 2 5- Right: Width and 
height are given by r c . The classifications are: NGC1530, TW; NGC4314, LW; NGC6300, 
CS; NGC2964, C. The white circles in both panels have radius r c . Dark regions are due to 
dust, while bright regions are due to emission. North is up and East is to the left. 
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Fig. 8. — Comparisons of barstrengths Qb and structure map rms a sm for SGD and LGD 
galaxies. The vertical widths correspond to the central 68-percentile spread of the cumulative 
fraction and are due to the uncertainty in Qb and <r sm . SGD galaxies are less strongly barred 
(left) and have less dust structure (right) than LGD galaxies. There is a 0.4% probability 
that the Qb are drawn from the same parent population, and 3% probability that the rms 
values are drawn from the same parent distribution. 
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Visually, SGD galaxies seem to have less circumnuclear dust than their LGD counter- 
parts. A K-S test comparing the <7 sm values of the two samples supports this impression, 
yielding a 3% probability that the two are drawn from the same parent distribution (see 
Figure 8). This is in agreement with the overall trend discussed below in §4.2 that more 
strongly barred galaxies have dustier circumnuclear regions, suggesting that the gas and dust 
funnelled to the center of a barred galaxy may eventually suffice to disintegrate a symmetric 
two-armed spiral in the central few hundred parsecs. The most likely mechanism for erasing 
the small-scale spiral structure is self-gravity, and later, energy feedback from substantial 
star formation (Wada et al. 2002, e.g.,). We note that many of the galaxies with the largest 
a sm (e.g., NGC1808, NGC2964, NGC6217 and NGC7552) are well-known starbursts or nu- 
clear starbursts. This indicates that there may be sufficient feedback from star formation to 
shape the structure of the circumnuclear ISM. 



4.2. Circumnuclear Rings and Strongly Barred Galaxies 

The material funnelled inward by the large-scale bar may not reach the central region at 
all, but instead form a circumnuclear starburst ring. If bars prompt rings of circumnuclear 
star formation, then star-forming rings should be more prevalent in strongly barred galaxies 
than in weakly barred ones. Observations do suggest that barred galaxies have a higher 
incidence of nuclear starburst rings than unbarred galaxies (see Kormendy & Kennicutt 
2004, and references therein). We offer here more evidence that this is in fact the case. 

Of the 75 galaxies studied, eight have circumnuclear starburst rings visible at radii 
less than or equal to 5% of D 2 5- All eight exhibit LGD structure, 3 as shown in Figure 10. 
One interesting transition case is NGC5427, where the LGD spiral arms appear to become 
the incomplete starburst ring itself. NGC4314 also exhibits this morphology, although the 
ring is more complete. The LGD dust lanes enter a number of other rings (e.g., NGC1097 
and NGC6951), but these "breaks" may simply be due to dust attenuation of a complete 
ring. As seen in Figure 1, the effective bulge radius and the circumnuclear ring radius for 
each of these galaxies apparently coincide. Of these eight galaxies, six have loosely wound 
circumnuclear spirals within the ring, while the other two (NGC864 and NGC1326) have 
chaotic circumnuclear spirals. Only 32 of the 75 galaxies are LW or CS; according to a 
binomial distribution, the probability of randomly selecting 8 galaxies that are restricted to 



3 The LGD structure in NGC1300 and NGC0864, while visible on the structure maps, is more prominent 
on the V — H color maps of Martini ct al. (2003a). This is because in these two cases, the relevant dust 
lanes are larger than the PSF, and thus have low contrast on the structure map. 
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these two classes is 0.062%. The circumnuclear morphology within a circumnuclear ring is 
therefore distinct from the ring itself. Visually, the rings appear highly turbulent and almost 
flocculent, while the dust structure in their interiors is much smoother. 

The galaxies with circumnuclear rings are also, as a group, more strongly barred than 
typical galaxies. The K-S test gives a probability of 1.5% that the ringed galaxies have the 
same barstrength distribution as the other 67 galaxies, and a probability of 3% that the LW 
ringed galaxies have the same distribution of barstrengths as the other fourteen LW galaxies. 
The barstrengths of the entire, LW, and ringed galaxy samples are shown in Figure 9. 
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Fig. 9. — Comparison of Qb for LW galaxies with and without circumnuclear rings. The 
vertical widths correspond to the central 68-percentile spread of the cumulative fraction and 
are due to the uncertainty in Q^. Galaxies with circumnuclear rings are found to be more 
strongly barred than typical galaxies, with 1.5% probability of being drawn from the same 
parent population. 
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There remain many pieces needed to assemble the circumnuclear starburst ring puzzle. 
A large-scale bar drives dust and gas towards the center of the galaxy. At the boundary 
of the bulge and the disk, a starforming ring develops. Within this ring, there is a fairly 
coherent loosely wound dust spiral. It is unclear, however, whether this nuclear dust spiral is 
composed of dust and gas originating from the large-scale bar and/or the ring itself, or if it is 
native dust that is "stirred-up" by the starburst ring. Simulations by Regan & Teuben (2004) 
and Maciejewski (2004) imply that bars in low velocity-dispersion systems are only capable 
of driving gas down to the radius at which the ring forms — but not further — suggesting that 
native dust comprises the nuclear spiral. A galaxy with a high central velocity-dispersion, 
however, could potentially have mass transported into the ring (Maciejewski 2004). 
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Fig. 10. — Structure maps of the eight galaxies in our sample with circumnuclear rings in 
order of increasing Q^. All eight have LGD structure. Width and height are given by O.I-D25. 
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For the entire sample of 75 galaxies, we find an increase in the structure map rms <7 sm 
with the barstrength Q&, although with substantial scatter and evidence the relation appears 
to turn over at large Qi>. The observed increase in dust structure may be due to increased 
star formation rates, which are correlated with the gas surface density via the global Schmidt 
law (Schmidt 1959; Kennicutt 1998). The structure maps for the twelve most strongly barred 
galaxies in our sample (Qb > 0.4) are shown in Figure 11. Of these twelve, five have LGD 
spirals, two of which end at a circumnuclear ring. The other seven have either chaotic or 
chaotic spiral nuclear dust morphology (5 C, 2 CS). As a circumnuclear ring is associated with 
nuclear star formation, and the chaotic structures seen in these strongly barred galaxies could 
potentially be related to star formation, it would be interesting to investigate the relationship 
between nuclear star formation rates, circumnuclear dust morphology, dust structure <j sm , 
and barstrength. 
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Fig. 11. — Structure maps for the 12 most strongly barred galaxies in our sample (Qb > 0.4). 
Qb increases to the right and then down. These galaxies have LGD structure (NGC1300, 
NGC1365, NGC5643, NGC4314, and NGC1530), chaotic nuclear structure (NGC6221, 
NGC3359, NGC1073, NGC4027, and NGC7479), or a chaotic nuclear spiral (NGC0613 and 
NGC2442). Two of the LGD spirals are associated with a circumnuclear ring (NGC1300 
and NGC4314). Dark regions are due to dust; bright regions are due to emission. 
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4.3. Large-scale Morphology: SB Galaxies 

One way in which barred spiral galaxies are commonly divided into subclasses is by 
taking note of where the large scale spiral starts with relation to the bar (e.g., Sandage 
& Bedke 1994). In SB(s) galaxies, the spiral arms begin at the ends of the bar, whereas 
in SB(r) galaxies, the spiral arms begin on a ring connecting the ends of the bar. SB(rs) 
galaxies are a transition group. SB(s) structures are thought to be preferentially found in less 
strongly barred galaxies than their SB(r) counterparts (e.g., Sanders & Tubbs 1980; Simkin 
et al. 1980). Furthermore, SB(s) galaxies typically show large-scale dust lanes that are not 
present in SB(r) galaxies, and SB(r) galaxies are observed to have less dust in their central 
regions than SB(s) galaxies (Kormendy & Kennicutt 2004). This is an apparent inconsistency 
with the generic bar-fueling picture: SB(r) galaxies are thought to have less dust — but be 
more strongly barred — while more strongly barred galaxies should have a higher central dust 
content. As these conclusions have been based on measuring the barstrength as an axis ratio 
(e.g., Sanders & Tubbs 1980), rather than as a force ratio, this discrepancy might be because 
SB(r) galaxies are actually more weakly barred than SB(s) galaxies. We reinvestigate the 
relation between these bar sub-types and barstrength. 

The breakdown by nuclear class, including LGD, of the 21 RC3-classified SB galaxies 
in our sample is given in Table 2. A small increase in Q b is seen from SB(r) to SB(s). 
Visually, the SB(r) galaxies have much less dust structure than the SB(rs) and SB(s) galaxies. 
According to the Wilcoxon test, the SB(r) sample has a smaller Qb than the SB(s) galaxies 
at a confidence level of 94%. Sanders & Tubbs (1980) suggest bar pattern speed as an 
alternative origin to the differences in SB(r) and SB(s) structure: a slowly rotating bar 
should give rise to SB(r) structure, while a rapidly rotating bar should yield SB(s) structure. 
As the differences in barstrength are reversed from what was expected, rotation may be more 
important in determining the large-scale morphology. With respect to the amount of dust 
structure, we find with > 99% confidence that the SB(r) galaxies have less dust structure 
(smaller a sm ) than the SB(s) sample; similarly, we find with 97% confidence that the SB(r) 
sample has less dust than the SB(rs) galaxies. 
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SB 


SB(r) 


SB(rs) 


SB(s) 


(Q b ) 


0.397 


0.272 


0.398 


0.450 


Total 


21 


3 


11 


7 


GD 


1 


1 








TW 


1 





1 





LW 


7 





5 


2 


cs 


3 





2 


1 


c 


8 


1 


3 


4 


N 


1 


1 








LGD 


6 





6 


2 



Table 2: Average barstrengths and nuclear classifications for SB galaxies. Nuclear classes 
are denned in §3.2, and LGD morphology is denned in §4.1. SB classification is given by the 
RC3 catalog. 
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sm 



Fig. 12. — Comparisons of structure map rms <r sm of SB galaxies. The vertical spreads 
correspond to the central 68-percentile spread of the cumulative fraction and are due to the 
uncertainty in Qt,. SB(r) galaxies have less dust structure than SB(rs) galaxies, which, in 
turn, have less dust than SB(s) galaxies. Specifically, with > 99% confidence, SB(r) galaxies 
have less central dust structure than SB(s) galaxies. 
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An important caveat for this aspect of our work is that large-scale morphological clas- 
sifications are known to depend on wavelength. For example, while only about one third of 
galaxies are classified as SB at visible wavelengths, over two thirds are classified as SB in 
the near infrared (Eskridge et al. 2000). This is primarily because near infrared light traces 
the underlying stellar mass, which is also one of the main reasons Q b is calculated using 
near infrared data (Buta & Block 2001). Eskridge et al. (2002) reclassify 205 spiral galaxies 
according to their iJ-band photometry from the OSUBGS. Forty-eight of our 75 galaxies 
are in this sample; 58% of which (28 of the 48) are classified as SB. Eskridge et al. do not, 
however, subclassify the if -band images into SB(r), (rs), and (s). It would be interesting to 
see how the dust content and Q b vary by class in NIR-identified SB galaxies. 



5. Conclusion 

We present a study of the circumnuclear dust morphology for a sample of 75 galaxies 
with archival HST data and measured barstrength Q b . This Q b is a measure of the maximal 
force ratio due to the presence of a bar in a galaxy, and thus is arguably superior to the 
rudimentary bar axis ratio (Buta & Block 2001). We use the structure map technique of 
Pogge & Martini (2002) to enhance the visibility of the galaxies' dust content and to classify 
the circumnuclear dust morphology within the central 5% of _D 25 , according to a refined 
version of the classification scheme proposed by Martini et al. (2003a). We also introduce the 
structure map rms <7 sm within the central regions of the galaxy as a quantitative measure of 
the amount of nuclear dust structure. A comparison of the morphological classifications and 
measured barstrengths reveals that tightly wound nuclear dust spirals (with pitch angles less 
than 10°) are preferentially found in galaxies with lower Q b . No other nuclear class exhibits 
a significant correlation with barstrength. 

Previous observations found that grand design nuclear dust spirals are hosted exclusively 
by strongly barred galaxies. While we do see grand design structure in many barred galaxies, 
this grand design structure does not extend all the way into the unresolved nucleus (~ 10 pc), 
although it is often present at larger scales. Earlier studies did not strictly require that the 
grand design structure extend into the nucleus. Taking this into account, we identify two 
distinct types of circumnuclear grand design spirals. Small grand design (SGD) spirals are 
nuclear dust spirals in which the two symmetric spiral arms are coherent from 1% of D 2 s 
(typically a few hundred parsecs) to the central tens of parsecs of the galaxy (where tracing 
the structure becomes resolution-limited). Large grand design (LGD) nuclear spirals, on the 
other hand, show two prominent symmetric arms within 10% of D25 (typically on the scale 
of a few kiloparsecs) ; these arms do not necessarily extend to the center of the galaxy. In 
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fact, these two types of grand design structure are nearly disjoint: the nuclear spiral arms 
in only two of the twenty LGD galaxies in our sample extend into the unresolved center of 
the galaxy (i.e., also displayed SGD structure). The LGD spirals are found in systematically 
more strongly barred host galaxies. This strongly confirms previous indications from much 
smaller samples and demonstrates that the dust lanes along the leading edges of large-scale 
spirals do not generally extend all the way into the nuclear region, but instead lose coherence 
at the scale of several hundred parsecs. While SGD spirals are not found in galaxies with 
a significantly different barstrength than typical galaxies, they are found in galaxies with 
significantly less dust structure in the central regions than LGD galaxies. The reduced 
dust structure may reflect a requirement for the formation of SGD morphology, or simply a 
requirement for its detection. 

The LGD spiral arms may not maintain coherence to the nucleus because mass inflow 
due to the presence of a large-scale bar prompts star formation, which can disrupt a nuclear 
grand design spiral. In addition, forty percent of the LGD spirals do not extend into the 
nuclear region because there is a circumnuclear starburst ring. We find that all of the galaxies 
with circumnuclear starburst rings have LGD structure and are more strongly barred than 
other galaxies. Within the rings, three fourths of the galaxies with circumnuclear rings have 
coherent loosely wound spirals within the ring, while the others have less coherent chaotic 
spirals. 

We also find that among SB galaxies, SB(s) galaxies have more dust structure and are 
more strongly barred than SB(r) galaxies. This is partially at odds with the prevailing view 
in the literature, which is that SB(r) galaxies should be more strongly barred, although 
there is consensus that SB(r) galaxies have less central dust (e.g. Kormendy & Kennicutt 
2004). Sanders & Tubbs (1980) suggest that differences in the bar pattern speed may also 
explain the large-scale morphological differences between SB(s) and SB(r) galaxies; as our 
results indicate that SB(s) galaxies are more strongly barred, pattern speed may be the more 
relevant parameter. 

Overall, there is agreement in the literature that more strongly barred galaxies should— 
and do — have more dust and gas in their centers (Kormendy & Kennicutt 2004). We find 
that for the most strongly barred galaxies, there are several possible morphologies the cir- 
cumnuclear dust can take. This may indicate that not all bars of a given strength Qb funnel 
material toward the centers of galaxies with equal efficiency, potentially due to the effects 
of pattern speed on bar efficiency, or simply the fact that Qi> is a one-parameter description 
of the bar. In the most strongly barred galaxies, there can be an LGD spiral whose arms 
do not extend to the galaxy nucleus but instead lose coherence. Some LGD spirals end 
in circumnuclear starburst rings. In the absence of these structures, however, the nuclear 
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dust in the most strongly barred galaxies tends to be fairly chaotic, potentially hosting star 
formation. It would be interesting to investigate how circumnuclear dust morphology and 
dust structure <r sm varies with nuclear star formation rates in strongly barred galaxies. 
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